Phonon dispersion curves and thermodynamic properties of C1-PU2O3 
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A recent inelastic x-ray scattering study [Manley et al., Phys. Rev. B 85, 132301 (2012)] reveals 
that the phonon dispersion curves of Pu02 is considerably consistent with our previous density 
functional +U results [Zhang et al., Phys. Rev. B 82, 144110 (2010)]. Here in the present work, 
using the same computational methods, we further obtain the phonon dispersion curves for Q-PU2O3. 
We find that the Pu-0 bonding is weaker in a-Pu203 than in fluorite PuC>2, and subsequently a 
frequency gap appears between the vibrations of oxygen and plutonium atoms. Based on the phonon 
dispersion curves and Helmholtz free energies of PuC>2 and a-Pu2 03, we systematically calculate the 
reaction energies for the transformations between Pu, PUO2, and a-Pu203. It is revealed that the 
thermodynamic equilibrium of the system is dependent on temperature as well as on the chemical 
environment. High temperature and insufficient oxygen environment are in favor of the formation 

Of CH-PU2O3. 

PACS numbers: 63.20.D-, 65.40.-b, 63.20.dk. 
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I. INTRODUCTION 

Plutonium-based materials attract much interest not 
only owing to their industrial, military, and environ- 
mental importance, but also for their basic theoretical 
prospects. Due to the complex character of Pu 5/ elec- 
trons, which locates in the boundary of localized and 
delocalized among the actinide metals, Pu occurs in six 
phases under different temperatures and pressures [l|, [|J . 
When metallic Pu is exposed to dry air at room tem- 
perature, the protective plutonium dioxide Pu02 layer is 
formed, which can further reduce to another thin layer 
of plutonium sesquioxide a-Pu2C>3 at the oxide- metal in- 
terface @- Since plutonium oxides can be applied to 
the nuclear stockpile and storage of surplus plutonium, 
revealing the physical properties of different plutonium 
oxides and the thermodynamic equilibrium criteria be- 
tween them become necessary and important. In partic- 
ular, the temperature dependence of the free energy of 
different plutonium oxides are key factors for their ther- 
modynamic equilibrium. 

The phase diagram of stoichiometric plutonium- 
oxygen system shows the presence of plutonium dioxide 
PuC>2 and sesquioxide PU2O3. PuC>2 shares the crys- 
tal structure of fluorite type, isomorphous with An02 
(An=\J, Th and Np). The actinide sesquioxide crystal- 
lizes into three different crystal structures, the hexago- 
nal La203 structure, the monoclinic Sn^Oa structure, 
and the cubic Mn 2 03 structure, respectively Q. PU2O3 
has been synthesized only in the Mn203 (a-Pu2 03) and 
La 2 3 structures (/3-Pu 2 3 ) Q. The a-Pu203 is sta- 
ble only below 300 °C. Structurally, the a-Pu203 unit 
cell can be obtained by ordered removal of 25 percent of 
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the oxide ions from PuC>2 2x2x2 supercell. Due to the 
structural similarity with PuC>2, a-Pu2 03 is more likely 
to be the reduction product of PuC>2. In fact, a mixture 
of these two oxides can be prepared by partial reduction 
of PUO2 at high temperature and then cooling to room 
temperature [5J. 

Because of its structural complexity (64 atoms in each 
unit cell), a-Pu2C>3 has seldom been theoretically stud- 
ied. Nevertheless, it is key to understand the thermo- 
dynamic equilibrium of the plutonium surface under dif- 
ferent oxidizing environments. Therefore in the present 
work, we systematically study the physical and thermo- 
dynamic properties of a-Pu203, the temperature depen- 
dent reaction energy between a-Pu203 and Pu02 is also 
discussed. 

Apart from the structural complexity of PU2O3, mod- 
eling of the electron localization/delocalization of pluto- 
nium 5/ electrons is also a complex task. Conventional 
density functional theory (DFT) schemes that apply the 
local density approximation (LDA) or the generalized 
gradient approximation (GGA) underestimate the strong 
on-site Coulomb repulsion of the plutonium 5/ electrons 
and consequently fail to capture the correlation-driven 
localization. However, the DFT+C/ method proposed by 
Dudarev et al. [1, 0] is an effective way to deal with 
the strong on-site Coulomb repulsion of Pu, U and Np 5/ 
electrons. In our previous studies, we have verified the 
validity and reliability of this method more than once [H- 
[ToL H2I [l3j , and the obtained phonon dispersion curves of 
Pu02 are very comparable with latest experiments [ill ]. 



II. CALCULATION METHODS 

The density functional theory (DFT) calculations are 
carried out usin g th e Vienna ah initio simulations pack- 
age (VASP) [l4 1 15l| with the projected-augmented-wave 
(PAW) potential method [lq |. The plane- wave basis 
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FIG. 1: (a) Atomic structure and magnetic arrangement of 
a-Pu2 03 within one unit cell, (b) Atomic structure of the 
fiuorite PuC>2. 

set is limited with an energy cutoff of 500 eV. The ex- 
change and correlation effects are described by gener- 
alized gradient approximation (GGA) in the Perdew- 
Burke-Ernzerhof (PBE) form [17| . The computational 
cells for a-Pu 2 03 and Pu0 2 contain 64 and 12 atoms re- 
spectively. And integration over their Brillouin Zones are 
done on 5 x 5 x 5 and 15x15x15 fc-ponit meshes generated 
by using the Monkhorst-Pack method [l8j], respectively. 
These fc-point meshes are carefully tested and proven to 
be sufficient for energy convergence of less than 1.0 xlO -4 
eV per atom. 

The strong on-site Coulomb repulsion among the lo- 
calized Pu 5/ electrons is described by the GGA+Z7 for- 
malisms formulated by Dudarev et al.$, 0|. As con- 
cluded in some previous studies, although the pure GGA 
or LDA fail to depict the electronic structure, especially 
the insulating nature and the occupied-state character of 
Pu0 2 and a-Pu 2 3 [1 [M El El, the LDA/GGA+t/ 
approaches will capture the Mott insulating properties 
of the strongly correlated Pu 5/ electrons in plutonium 
oxides adequately. Since only the difference between the 
spherically averaged screened Coulomb energy U and the 
exchange energy J is significant for the total LDA/GGA 
energy functional, we labeled them as one single effective 
parameter U for simplicity. The value of exchange energy 
J is set to be a constant 0.75 eV, which is the same as in 
our previous study of plutonium oxides, [il. ITol. [l3j while 
U can suitably be chosen within a range of values (~1- 
4 eV), which are in the ballpark of commonly accepted 
values for actinide elements. 



III. RESULTS AND DISCUSSIONS 

A. Ground-state properties of Q-PU2O3 and P11O2 

For GGA/LDA, Coulomb repulsion of the 5/ electrons 
in Pu describes Pu0 2 as incorrect ferromagnetic (FM) 
conductor instead of antiferromagnetic (AFM) Mott in- 
sulator reported by experiment [20J . Our previous studies 



show that the DFT+J7 approach can accurately describe 
the atomic and electronic structures of plutonium dioxide 
d, EH . The AFM phase begins to be energetically pre- 
ferred at Hubbard J7~1.5 eV, consequently, AFM phase 
becomes more stable than FM phase above U =2 eV. As 
shown in Fig. 1, we take use of the collinear AFM or- 
der for a-Pu20,3 unit cell and 2x2x2 Pu0 2 supercell 
as proposed by Regulski et al. 21], since we have con- 
firmed that this spin order is energetically favored fl3j . 
We reported in Table I our calculated results of lattice 
parameter ao, bulk modulus B<j, pressure derivative of 
the bulk modulus B' , spin moments /i mas , and energy 
band gap E g for AFM a-Pu 2 3 and Pu0 2 at GPa 
within GGA and GGA+L/ (t/=3 eV). For comparison, 
experimental values in Refs. 0, H, HH, [23[ are also listed. 
All the values of ao, -Bo, and B' Q are obtained by fitting 
the third-order Brich-Murnaghan equation of state [24j. 
The lattice parameter ao of a-Pu 2 03 is 11.169 A at U=3 
eV, in accordance with the measured value range 11.03- 
11.07 A [I, H|. Comparing the volumes of Pu0 2 2x2x2 
supercell and a-Pu 2 03 unit cell, there is about 7% ex- 
pansion when the reduction process accomplishes. The 
bulk modulus Bo of a-Pu 2 03 are calculated to be 120 
GPa, which is much lower than the 192 GPa value of 
Pu0 2 , indicating that the resistance to external pressure 
of a-Pu 2 03 is weaker than Pu0 2 . 

An accurate description of the electronic structure for 
materials is of great importance. The photoemission ex- 
periments show that the Pu0 2 is Mott insulator [25|, [26[ . 
In our previous studies, the insulating nature and the 
occupied-state character of plutonium dioxide, which the 
pure GGA fail to depict, can be prominently improved 
by turning the effective Hubbard U parameter in a rea- 
sonable range. When we turned on U parameter, the 
amplitude of insulating gap increases with enhancing U 
for a-Pu 2 03. Similar to Pu0 2 , the FM conductor to 
AFM insulator transition of a-Pu 2 03 also occurs, and 
at U—3 eV the insulating gap is 1.0 eV. The calculated 
amplitude of local spin moment is 5.02 /is per Pu atom 
for AFM a-Pu 2 03 within GGA+U scheme, somewhat 
larger than 4.14 [Ib of Pu0 2 . 



B. Phonon dispersions, thermal expansion, and 
heat capacity of Q-PU2O3 

The phonon dispersion curves of a material is closely 
related to its structure stability and thermodynamic 
properties like thermal expansion, Helmholtz free energy, 
and heat capacity. In our present study, the phonon fre- 
quencies are carried out using the direct method. For cal- 
culating the Hessian matrix of a-Pu 2 03 we adopted its 
unit cell with the size of over 11 Aalong the x,y,z direc- 
tions, which is large enough to isolate direct interactions 
of neighboring cells. And a 3x3x3 Monkhorst-Pack k- 
point mesh is used for integration over the Brillouin Zone 
while calculating the forces. The calculated phonon dis- 
persion curves of a-Pu 2 03 along the high-symmetry lines 
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TABLE I: Calculated lattice parameters ao, bulk modulus Bo, 
pressure derivative of the bulk modulus B' , spin moments 
fJ-mag, and energy band gap E g for AFM a-Pu2C>3 and PuC>2 
in the GGA and GGA+U ((7=3 eV) approaches at GPa. 
For comparison, experimental values are also listed. 



Method 


ao 


Bo 


B' 


[l>rnag 


E 9 




(A) 


(GPa) 




(Ms) 


(eV) 


a-Pu 2 3 GGA 


10.921 


123 


4.25 


4.95 


0.0 


GGA+U 


11.169 


120 


3.53 


5.02 


1.0 


Expt. 


11.03-11. 07 a - b 










PuC-2 GGA 


5.396 


185 


4.29 


4.09 


0.0 


GGA+U 


5.457 


192 


4.50 


4.14 


1.2 


Expt. 


5.398 c 


178 d 









a Reference 
b Reference 
c Reference 
d Reference 
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FIG. 4: Temperature dependence of the bulk modulus B for 
a-Pu2 03. The inset is the ratio of B/Bo for a-Pu2 03 and 
PuQ 2 . 
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FIG. 5: Heat capacities of a-Pu203 and Pu02 at constant 
volume (CV) and constant pressure (C p ). 



FIG. 2: Phonon dispersion curves (left) and phonon density 
of states (Phonon DOS) (right) of a-Pu 2 3 . 
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FIG. 3: Temperature dependence of the linear thermal expan- 
sion for a-Pu203 and Pu0 2 . The inset is the volume thermal 
expansion coefficient as a function of temperature. 



of its Brillouin Zone are shown in Fig. 2, together with 
the atomic phonon density of state (phonon-DOS). As 
presented, there is an evident frequency gap between the 
oxygen and plutonium vibrations. The vibrations of the 
plutonium and oxygen atoms dominate the low (0-5.2 
THz) and high frequency (8-16.5 THz) parts respectively. 
The obvious frequency gap between oxygen and pluto- 
nium vibrations makes a-Pu2 03 different from PuC>2, as 
there is a small overlap between oxygen and plutonium 
vibration modes in Pu0 2 [HI E3 • 

From Fig. 2 we see that the frequency value of the 
phonon dispersion curves for a-Pu2C>3 can only be as 
large as 16.63 THz, which is much smaller than the max- 
imum value of 19.19 THz for PuC>2. Since the high- 
frequency part corresponds to oxygen vibrations, this re- 
sult indicates that the binding of oxygen atoms around 
plutonium atoms is stronger in Pu0 2 than in Q-PU2O3. 
In addition, as mentioned above, the C1-PU2O3 unit cell 
can be obtained by orderly removing 25 percent of the 
oxide atoms from a 2x2x2 PuC>2 supercell. The re- 
moval of oxygen atoms leads to a 7% volume expansion 
of a-Pu2 03, and then weakens the bond strength of Pu- 
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O bondings. Overall, the all-positive phonon dispersion 
curves confirm the thermodynamic stability of a-Pu203 
at its ground state. 

Based on the obtained phonon dispersion curvs, we 
can further calculate the thermodynamic properties of 
C11-PU2O3, which begins by calculating its Hclmholtz free 
energy F(V,T). Normally the of Helmholtz free energy 
of an insulator can be divided into two parts as 



F(V,T) = E(V) + F vib (V,T), 



(1) 



where E(V) is the ground-state electronic energy, and 
F vib (V, T) stands for the phonon free energy at a given 
temperature T. Within quasi-hamonic approximation 
(QHA), the F V ib(V,T) term can be evaluated by 



F vib (V,T) = k B Tj2^ 



2 sinh 



(2) 



where u)j (q, V) is the phonon frequency of the j th phonon 
mode with wave vector q at fixed cell volume V, and k B 
is the Boltzmann constant. The total specific heat of the 
crystal is the sum of all phonon modes over the BZ, 



C v {T) = Y,cvA^T). 



(3) 



cv,j(q,T) is the mode contribution to the specific heat 
defined as, 



c yj (q,T) = k B Y, 



fa; 3 -(q,vo y 1 

2k B T J sinh 2 [^(q,F)/2fc B r]' 

(4) 

From the obtained phonon dispersion curves and Eqs. 
(1) and (2), the Helmholtz free energy of a-Pu2C>3 is 
calculated at different temperatures and different lattice 
constants. By determining the minimum energy lattice 
constants at different temperatures, the lattice expansion 
curve is obtained and shown in Fig. 3. The correspond- 
ing result of PuC>2 is also shown in Fig. 3 for comparisons. 
One can see that the rate of lattice expansion is larger for 
a-Pu 2 3 than for Pu0 2 , indicating that the volume of a- 
PU2O3 is more sensitive to temperature. This character 
can also be seen from the thermal expansion coefficient 
ay (T) as shown in the inset of Fig. 3. The bulk modulus 
B is also calculated at different temperatures according 
to the formula 

d 2 F 

B = ^0(^2)1^=^0- 

For a-Pu203, the bulk modulus decreases with increas- 
ing temperature, and the decreasing rate is much larger 
than that of Pu02, which can be seen from the ratio of 
the bulk modulus B/Bo as depicted in Fig. 4. Within 
the framework of QHA, the considered vibration modes 
are harmonic but volume dependent. This volume de- 
pendence is somehow neglected in the above calculations 
because the considered lattice expansions are very small. 



By using Eqs. (3) and (4), the heat capacity at con- 
stant volume Cy can be calculated for CH-PU2O3 and 
Pu02, which are shown in Fig. 5. Their heat capac- 
ity at constant pressure C p is also calculated from the 
relationship 

C p - C v = a 2 v (T)B(T)V(T)T 

, and shown in Fig. 5. We can see that as temperature 
increases, the values of Cy and C p increase continuously 
with the growth rate of Cy slightly less than C p , both 
for a-Pu2 3 and Pu02- Comparing the specific heat 
of a-Pu2 3 and Pu02, there is a crossing point at 265 
K. Beneath 265 K, the increasing rate of heat capacities 
for a-Pu203 is larger than that for Pu02- At the room 
temperature of 300 K, the values of Cy and C p for a- 
Pu 2 3 become 196.1 J kg^ 1 K" 1 and 199.0 J kg" 1 K" 1 , 
respectively. 



C. Thermodynamic equilibrium of different 
plutonium oxides 

Understanding the thermodynamic equilibrium of dif- 
ferent plutonium oxides at different temperatures can 
help people get more clear about the oxidation of plu- 
tonium, and further enhance nuclear security and ma- 
nipulate the safe maintenance of plutonium. Chemically, 
the reaction energy AF can be defined as the energy 
differences between initial (reactant) and final (product) 
states of a reaction, i.e., AF = F product ~ F reactant . If 
a reaction releases energy (AF<0), it is thcrmodynami- 
cally favorable and can occur spontaneously. Conversely, 
if AF is positive, the reaction is not spontaneous and is 
hindered until sufficient energy is added to the system. 

By interacting with oxygen molecules, plutonium can 
be oxidized according to the processes 



and 



Pu + 2 ^Pu0 2 , 



Pu+^0 2 -> ipu 2 3 . 



(5) 



(6) 



We have known that both a-Pu203 and Pu02 are stable 
in solid phases at room temperature. As the physical 
and chemical conditions of the plutonium surface change, 
these two oxides are able to convert to each other. When 
excess O2 is present, a-Pu203 is further oxidized through 
the process 



Pu 2 3 + io 2 -> 2Pu0 2 . 



(7) 



In contrast, if oxygen is sufficient in the environment, 
Pu02 can be reduced to a-Pu2 03 through the process, 



3Pu0 2 +Pu -> 2Pu 2 3 . 



(8) 



Based on their structural similarities, here we only fo- 
cus on the transformation between Pu02 and CH-PU2O3 
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TABLE II: Reaction energies AF for the formation of different plutonium oxides at different temperatures. The parenthesized 
(Pu) and (O2) represent for the reaction energy for 1 mole Pu atoms and 1 mole O2 molecules, respectively. For comparison, 
the experimental data in Ref. are also listed. 





K 




300 K 




500 K 




Expt.* 




Reaction 


kcal/mol 


kcal/mol 


kcal/mol 


kcal/mol 


kcal/mol 


kcal/mol 


kcal/mol 


kcal/mol 




(Pu) 


(O2) 


(Pu) 


(O2) 


(Pu) 


(O2) 


(Pu) 


(O2) 


Pu+0 2 -> Pu0 2 


-211.0 


-211.0 


-213.4 


-213.4 


-215.5 


-215.5 


-239 


-239 


Pu+f 2 -> iPu 2 3 


-176.0 


-234.7 


-178.0 


-237.3 


-181.5 


-242.0 


-189 


-252 


Pu 2 3 +i0 2 2Pu0 2 




-139.9 




-137.8 




-136.1 






3Pu0 2 +Pu 2Pu 2 3 


-71.1 




-74.6 




-79.4 









* Reference^ 



without considering the (3 phase of PU2O3. During sur- 
face oxidation of plutonium, temperature plays a signifi- 
cant role for the system to reach thermodynamic equilib- 
rium. For this reason, we calculated the reaction energies 
for all the above four processes at different temperatures, 
and the results are listed in Table II. As a comparison, 
the experimentally determined values by Haschke et al. 
[H are also included. In our calculations, bulk S-Pu is 
chosen to be the chemical source of Pu, and its energy 
per Pu atoms represents for the chemical potential of 
Pu. The chemical potential of oxygen is chosen to be 
the cohesive energy of an gas-like O2 molecule. To avoid 
the overestimation of DFT methods on the 0-0 bond- 
ing strength, we adopt the experimental cohesive energy 
of 5.21 eV [23| . From our calculated results, one can see 
that all the four reaction processes have negative reaction 
energies, indicating that they are thermodynamically al- 
lowable. We can also see from Table II that as temper- 
ature increases, the energy release for reactions (5), (6), 
and (8) increases, indicating that these three reactions 
can be accelerated by higher temperatures. Especially, 
we can see that formation of PU2O3 is more efficient at 
higher temperatures. In contrast, the oxidation reaction 
of PU2O3 (reaction (7)) releases more energy at lower 
temperatures. Thus fully oxidation of 5-Pu into Pu02 is 
more efficient at lower temperatures. 

The reaction energy values per mole Pu correspond 
to the situation where the quantity of Pu is fixed and 
oxygen is sufficient, while the values per mole O2 cor- 
respond to the situation where the quantity of oxygen 
is fixed and Pu is sufficient. We can see from Table II 
that the energy release per mole Pu for oxidizing <5-Pu 
into PUO2 is always larger than into a-Pu203, indicating 
that formation of Pu02 is more favorable when oxygen is 
sufficient in the environment. Correspondingly, the en- 
ergy release per mole O2 for oxidizing (5-Pu into Pu02 is 
always smaller than into a-Pu203, indicating that forma- 
tion of CV-PU2O3 is more favorable when Pu is sufficient 
in the environment. This conclusion accords well with 
the experimental results. 

IV. SUMMARY 

In summary, we have performed a systematical den- 
sity functional theory study to investigate the thermody- 



namic properties of a-Pu203 and Pu02, as well as the 
thermodynamic equilibrium between them. The antifer- 
romagnetic order and insulator character of a-Pu203 can 
be well predicted by the GGA+U (U=3 eV) approach. 
The theoretical atomic and electronic structures of the 
plutonium oxides are in agreement with available exper- 
imental results. Based on the calculated phonon dis- 
persion curves of a-Pu203 and Pu02, we systematically 
obtain their thermodynamic properties like thermal ex- 
pansion coefficient, bulk modulus, and heat capacities at 
different temperatures. We find that the volume expan- 
sion and bulk modulus of a-Pu203 are both more sensi- 
tive to temperature than Pu02, and there is a transition 
temperature of 265 K comparing the specific heats of a- 
P112O3 and Pu02- In the temperature range between 
OK and 265 K, a-Pu203 presents higher heat capacities 
than P11O2. The calculated reaction energies for the four 
processes during surface oxidation of 5-Pu suggest that 
the reaction products depend heavily on the tempera- 
ture and chemical condition of the plutonium surface. 
The theoretical reaction energies of Pu oxidizing into a- 
PU2O3 and PUO2 are in accord with experimental data 
and can reflect the basic qualitative characteristic of the 
real physical process. Besides, temperature plays a sig- 
nificant role in the component of oxidation product. Our 
results suggest that when temperature increases, Pu02 
tends to be reduced to be a-Pu2 03. 
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